Particle Acceleration in Turbulence and Weakly Stochastic Reconnection 
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Fast particles are accelerated in astrophysical environments by a variety of processes. Acceleration 
in reconnection sites has attracted the attention of researchers recently. In this letter we analyze 
the energy distribution evolution of test particles injected in three dimensional (3D) magnetohydro- 
dynamic (MHD) simulations of different magnetic reconnection configurations. When considering 
a single Sweet-Parker topology, the particles accelerate predominantly through a first-order Fermi 
process, as predicted in previous work 0] and demonstrated numerically in [l^- When turbulence 
is included within the current sheet, the acceleration rate, which depends on the reconnection rate, 
is highly enhanced. This is because reconnection in the presence of turbulence becomes fast and 
independent of resistivity [l^ . [Tgj and allows the formation of a thick volume filled with multiple 
simultaneously reconnecting magnetic fluxes. Charged particles trapped within this volume suffer 
several head-on scatterings with the contracting magnetic fluctuations, which significantly increase 
the acceleration rate and results in a first-order Fermi process. For comparison, we also tested 
acceleration in MIfD turbulence, where particles suffer collisions with approaching and receding 
magnetic irregularities, resulting in a reduced acceleration rate. We argue that the dominant accel- 
eration mechanism approaches a second order Fermi process in this case. 
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Introduction. Energetic particles are ubiquitous in as- 
trophysical environments. Cosmic ray (CR) acceleration 
still challenges the researchers. For instance, the origin 
of ultra high energy cosmic rays (UHECRs) is still un- 
known, although several mechanism have been proposed 
[29I and references there in]. Their spectrum is consis- 
tent with an origin in extragalactic astrophysical sources 
and candidates range from the birth of compact objects 
to explosions related to gamma-ray bursts (GRBs) or 
events in active galaxies (AGNs) Very high en- 

ergy observations of AGNs and GRBs with the Fermi and 
Swift satellites and ground based gamma ray observato- 
ries (HESS, VERITAS and MAGIC) challenge current 
theories of particle acceleration, mostly based on the ac- 
celeration in shocks, trying to explain how particles are 
accelerated to energies above TeV in regions relati vely 
small compared to the fiducial scale of their sources [33 ■ 

While particle acceleration in shocks has been exten- 
sively explored [s^] (see also 14 , 27 1 for reviews) , an alter- 
native, less investigated mechanism so far, involves parti- 
cle acceleration in magnetic reconnection sites. Magnetic 
reconnection may occur when two magnetic fluxes of op- 
posite polarity encounter each other. In the presence of 
finite magnetic diffusivity the converging magnetic lines 
annihilate at the discontinuity surface and a current sheet 
forms there. In 0] , the authors first proposed that an effi- 
cient first-order Fermi process can occur within a current 
sheet, where trapped charged particles may bounce back 
and forth several times and gain energy due to head- 



on collisions with the two converging magnetic fluxes in- 
coming with the reconnection speed Vrec- They found 
that the particle energy gain after each round trip is 
t^E jE oc Vrecjc. With a fast magnetic reconnection, like 
the one induced in the presence of turbulence [l9j . Vre^ 
can be of the order of Alfven speed Va- Afterwards, [3] 
appealed to a similar process, but within a collisionless 
reconnection scenario. In their model, the acceleration 
is controlled by the contraction of two-dimensional (2D) 
loops due to firehose instability that arises in a particle- 
in-cell (PIC) domain [see also lol. [26j . Other processes of 
acceleration, like those due to the electric field associ- 
ated with the current in the reconnection region |24l| and 
turbulence arising as a result of reconnection [l8|, were 
shown to be less dominant. 

Magnetic reconnection is expected to induce acceler- 
ation of particles in a wide range of galactic and ex- 
tragalactic environments. Originally discussed predomi- 

nantly in the context of solar fiares [e.g. [3-[i [li, [H [2O, 

21, 38], it has been recently extended to more extreme 
astrophysical environments and sources, such as the pro- 
duction of UHECRs Il4j . in particle acceleration in 
jet- accretion disk systems |4l46l. |ll|. and in the general 
framework of AGNs and GRBs [lllll [sl-lil] • These ap- 
plications, however, still require in-depth studies of par- 
ticle acceleration in magnetic reconnection sites, and its 
connection with magnetohydrodynamical (MHD) turbu- 
lence and fast magnetic reconnection. 

In this letter, we explore this issue by means of fully 
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3D MHD simulations of particle acceleration in magnetic 
reconnection domains. A p reliminary study in this direc- 
tion was performed in [l5| where the analytical model of 
Q was successfully tested. It was also shown that particle 
acceleration taking place in MHD reconnection domains 
without including kinetic effects produces results similar 
to those found in collisionless PIC simulations Q. This 
proved that the acceleration in reconnection regions is 
a universal process which is not determined by the de- 
tails of the plasma physics and can be very efficient in 
collisional gas, although energy and radiative losses due 
to the interactions of the accelerated particles with the 
surrounding plasma may be significant in some systems. 
Moreover, it has been shown that particle acceleration in 
2D and 3D MHD reconnection behaves quite differently, 
what calls for focusing on realistic 3D geometries. 

In this work we study differences between energy dis- 
tributions of test particles injected in 3D MHD reconnec- 
tion sites with and without turbulence, corresponding to 
the Sweet-Parker 



I34| and Lazarian-Vishniac [19[ con- 
figurations. For comparison, we also considered a pure 
turbulent environment. 

Numerical Simulations of Reconnection and Turbu- 
lence. Magnetic reconnection and turbulence were mod- 
eled by solving the MHD equations numerically on a 
uniform mesh using a shock-capturing Godunov-type 
scheme based on the 2"'' order spatial reconstruction 
and 2"*^ Runge-Kutta (RK) time integration [l^. We 
used the isothermal HLLD Riemann solver 28| to obtain 
the numerical fluxes in the update step, what resulted 
in reducing the energy dissipation of Alfven waves. We 
incorporated the field-interpolated constrained transport 
(CT) scheme based on a staggered mesh [2^ into the 
integration of the induction equation to maintain the 
V • B = constraint numerically. 

In the reconnection models we adopted an initial Har- 
ris current sheet, B^ix^y, z) = Boa; tanh(y/0), initialized 
by the magnetic vector potential with a uniform guide 
field Bz{x,y,z) = Bqz = const and the density profile 
set from the uniform total (thermal plus magnetic) pres- 
sure pT{x,y,z) = const. Initial velocity was zero every- 
where. The reconnection was initiated by a small initial 
perturbation 5Az{x,y,z) = (^Boa, cos(27ra;) exp[— (y/c?)^] 
added to the initial vector potential. The parameters 
SBqx and d describe the strength of the initial perturba- 
tion and thickness of the perturbed region, respectively. 
In the case of turbulent model the initial magnetic field 
was uniform. 

We employed dimensionless equations, so that the ve- 
locity and magnetic field are expressed in the fiducial 
Alfven speed units (defined by the antiparallel compo- 
nent of the field) and the unperturbed density po = !)• 
The length of the box in the X direction defines the dis- 
tance unit and the time is measured in units of L/Va- 
Initially, we set the strength of the antiparallel magnetic 
field component B^^ to 1.0 and the guide field Bqz to 



0.1. The sound speed was set to 4.0 in all models in or- 
der to suppress the compressibility in the system. The 
parameters describing the initial perturbation are set to 
5Bqx = 0.05 and d = 0.1. For more details of the com- 



putational method and numerical setup see [16|, |17| . 

Integration of Particle Trajectories. In order to in- 
tegrate the test particle trajectories we freeze in time a 
data cube obtained from the MHD models and inject test 
particles in the domain with random initial positions and 
directions and with an initial thermal distribution. For 
each particle we solve the relativistic motion equation 



— (jmu) = q(E + ux B) 
dt 



(1) 



where to, q and u are the particle mass, electric charge 
and velocity, respectively, E and B are the electric and 
magnetic fields, respectively, 7 = (l — /<?) is the 
Lorentz factor, and c is the speed of light. The electric 
field E is taken from the MHD simulations 
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-V X B + i]J, 



(2) 



where v is the plasma velocity, J = V x B is the current 
density, and 77 is the Ohmic resistivity coefficient. We ne- 
glect the resistive term above since its effect on particle 
acceleration is negligible Ts"]. In the current studies we 
do not include the particle energy losses, thus particles 
can gain or loose energy only through the interactions 
with the moving magnetized plasma. The inclusion of 
radiative losses or back reaction on the plasma is planned 
for future studies, particularly with the incorporation of 
the kinetic effects. We note that since we are focusing 
on the acceleration process only, we consider very sim- 
ple domains which represent only small periodic boxes 
of entire magnetic reconnection or turbulent sites. For 
this reason, the typical crossing time through the box of 
an injected thermal particle is very small and it has to 
re-enter the computational domain several times before 
gaining significant energy by multiple scatterings. Thus, 
whenever a particle reaches the box boundary it re-enters 
in the other side to continue scattering. 

Eq. [T] is integrated using the 4*'* order Runge-Kutta 
method with the adaptive time step based on the 5*'' 
order error estimator [see WH e.g.]. The fields are inter- 
polated using the 2"^^ order methods [2^ and limited to 
Qth Qj.(jgj. near discontinuities [is} . 

For convenience, we assume the speed of light to be 20 
Va, which defines our plasma in a non-relativistic regime, 
and the mean density is assumed to be 1 atomic mass 
unit per cubic centimeter, which is a fiducial value, e.g., 
of the interstellar medium (ISM) density. All times are 
expressed in units of the Alfven time. 

Acceleration in Sweet-Parker Reconnection. The ac- 
celeration of a single test particle in the Sweet-Parker 
reconnection was already described in [3]. Here, we 
present statistical studies for 10,000 protons injected in 
such a domain. 
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In the Sweet-Parker model [30|, |3J| the reconnection 
speed is given by Vrec ~ VaS~^^'^, where S = LVa/t] is 
the Lundquist number. Because of the typical huge as- 
trophysical sizes (L), S is also huge for Ohmic diffusivity 
values (e.g., for the ISM, S - lO^^). The Sweet-Parker 
reconnection is very slow, unless we use an artificially 
large diffusivity. In the model shown in the top panel of 
Figure [T] we employed a diffusivity coefficient 77 = 10"'^ 
expressed in code units. This value, due to the numeri- 
cal diffusivity, is several orders of magnitude larger than 
typical Ohmic diffusivity in astrophysical environments, 
and makes the Sweet-Parker reconnection in the simu- 
lation efficient. The time evolution of the energy dis- 
tribution for the accelerating particles is shown for this 
model in the top left panel of Figured] Initially, the per- 
pendicular acceleration dominates, because the volume 
in which we inject particles is much larger than the cur- 
rent sheet region (shown in the right panel of the same 
figure). The perpendicular acceleration, due to a drift of 
the magnetic flux, starts before the particles reach the 
reconnection region The distribution of particles 

does not change significantly until t = 1.0. Then, a rapid 
increase in energy by roughly four orders of magnitude 
appears for a fraction of particles. We observe a big gap 
between the energy levels before and after these acceler- 
ation events, which is also evident in the particle energy 
spectrum depicted in the embedded subplot of the same 
diagram. The events are spread in time because parti- 
cles gain substantial energy at different moments when 
crossing the current sheet. The energy growth during 
this stage is exponential. This is a first-order Fermi ac- 
celeration process, as predicted in [3|] and tested already 
in 
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Around t = 10, we see a transition from the exponen- 
tial to power-law particle energy dependence with esti- 
mated index a ^ 1.1 (see the upper plot of Figure [T]). 
A similar transition was also observed in [l5| . In the 
plot we show the kinetic energy normalized by the proton 
rest mass, what is equivalent to (7 — 1). The exponen- 
tial acceleration stops right after the energy value 10^ is 
reached, because the gyroradii exceed the thickness of the 
acceleration region. From this moment on the particles 
cannot be confined within this thickness and the first or- 
der Fermi process ceases. Further energy increase is due 
to a much slower drift acceleration (of the perpendicu- 
lar component only) caused by the large scale magnetic 
fields gradients. The presence of a guide field allows the 
particles to accelerate in the parallel direction as well. 

Although the Sweet-Parker model with an artificially 
enhanced resistivity results in a predominantly first-order 
Fermi acceleration, only a small fraction of the injected 
particles is trapped and efficiently accelerated in the cur- 
rent sheet (see the energy spectrum of the accelerated 
particles in the bottom right of Figure [IJ , because the 
acceleration zone is very thin. 



Acceleration in Reconnection with Turbulence. Lazar- 
ian & Vishniac proposed a model for fast reconnec- 
tion that does not depend on the magnetic diffusivity (see 
also [l3|)- The model appeals to the ubiquitous astro- 
physical turbulence as a universal trigger of fast recon- 
nection. The predictions of this model have been suc- 
cessfully tested in numerical simulations 1^ 12 1 which 
confirmed that the reconnection speed is of the order of 
the Alfven speed. An important consequence of the fast 
reconnection by turbulent magnetic fields is the forma- 
tion of a thick volume filled with small scale magnetic 
fluctuations. In order to test the acceleration of particles 
within such a domain, we introduced turbulence within 
a current sheet with a Sweet-Parker configuration (as de- 
scribed in the previous paragraph), then we followed the 
trajectories of 10,000 protons injected in this domain. 

The middle left panel of Figure [T] shows the evolution 
of the kinetic energy of the particles. After injection, a 
large fraction of test particles accelerates and the parti- 
cle energy growth occurs earlier than in the Sweet-Parker 
case (see also the energy spectrum at < = 5 in the detail 
at the bottom right of the same diagram). This is ex- 
plained by a combination of two effects: the presence of 
a large number of converging small scale current sheets 
and the broadening of the acceleration region due to the 
turbulence. Here, we do not observe a gap seen in the 
Sweet-Parker reconnection, because particles are contin- 
ually accelerated by encounters with several small and 
intermediate scale current sheets randomly distributed in 
the thick volume. The acceleration process is clearly still 
a first order Fermi process, as in the Sweet-Parker case, 
but more efflcient as it involves larger number of particles, 
since the size of the acceleration zone and the number of 
scatterers have been naturally increased by the presence 
of turbulence. Moreover, the reconnection speed, which 
in this case is independent of resistivity [ll,ll3 and de- 
termines the velocity at which the current sheets scat- 
ter particles, has been naturally increased as well (i.e. 

Vrec Va). 

During this stage a is in the range 2.48 — 2.75. Then, 
like in the laminar case, the protons accelerate at smaller 
rates after reaching energy level ^ W^, because the thick- 
ness of the acceleration region (and of the plasma scat- 
tering centers) becomes smaller than their Larmor radii. 
This, however, occurs later than in the previous case 
(around t = 100). 

Acceleration in Pure Turbulence. For comparison, in 
the bottom left panel of Figured] we show the kinetic en- 
ergy evolution of accelerated particles in a domain with 
turbulence injected in an initially uniform large scale 
field. In this domain the acceleration is initially less effi- 
cient and a much smaller fraction of particles is acceler- 
ated than in the reconnection case. The increase in the 
acceleration rate (with a ~ 2.54 — 3.21) is in part due to 
constraints of the computational domain. We cannot in- 
ject turbulence at scales larger than the size of the box. 
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FIG. 1. Left column: Particle kinetic energy distributions for 10,000 protons injected in the Sweet-Parker reconnection (top), 
fast magnetic reconnection (middle), and purely turbulent (bottom) domains. The colors indicate which velocity component 
is accelerated (red or blue for parallel or perpendicular, respectively). The energy is normalized by the rest proton mass. 
Subplots show the particle energy distributions at t = 5.0. Right column: The exemplary XY cuts through the domain at 
Z = of the absolute value of current density \J\ overlapped with the magnetic vectors for the Sweet-Parker reconnection 
(top), fast reconnection (middle), and purely turbulent domains (bottom). Models with Bqz = 0.1, ?7 = 10~^, and the resolution 
256x512x256 for reconnection and Boz = 0.2 and the resolution 128x256x128 for the turbulent cases are shown. 



As a consequence, an undesired converging flow arises 
along the mean field due to large scale Alfven waves, 
which enhances the acceleration rate for particles with 
Larmor radii approaching the turbulent injection scale. 
After reaching an energy level of about 10^ proton mass, 
this acceleration is significantly suppressed and a drops 
down to ^ 0.67. In the model with fast reconnection 
the presence of the large scale current sheet provides the 
converging flow. This flow brings scattering centers al- 
lowing a continuous growth of the particle energy until 
the saturation level. In pure turbulence, the absence of 
a converging flow results in a random particle scattering 
on approaching and receding small scale current sheets 
(although at a smaller rate), so that the overall acceler- 
ation is a second-order Fermi process. This point still 



requires further studies, as reconnection layers in pure 
turbulence can be responsible for first-order Fermi ac- 
celeration of low energy particles. As before, the rapid 
transition to smaller acceleration rate occurs when the 
particle gyroradius reaches the size of the turbulent do- 
main and its irregularities (around t — 10). We note that 
we have neglected here the time evolution of the MHD 
environment since this is much longer than the particle 
time scales. In fact, particles are accelerated by magnetic 
fiuctuations in the turbulent field and interact resonantly 
with larger and larger structures as their energy increases 
due to the scatterings. In a steady state turbulent envi- 
ronment, as considered here, particles will see on average 
the same sort of fluctuation distribution, so that after 
several Alfven times, one should expect no significant 
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changes in the particle spectrum due to the evolution 
of the large scale MHD environment. Nonetheless, this 
evolution may be important when considering more real- 
istic non-steady environments and when calculating real 
spectra and loss effects [e.g. [2^. 

Conclusions. In this letter we investigated particle 
acceleration in 3D MHD domains of magnetic recon- 
nection. We found that the presence of turbulence sig- 
nificantly increases the acceleration rate in a first-order 
Fermi process as predicted in Q. The particles trapped 
within the current sheet suffer several head-on scatterings 
with the contracting magnetic fluctuations in the thick 
volume. In the Sweet-Parker model, where the reconnec- 
tion speed was artificially large due to numerical mag- 
netic diffusivity, the acceleration rate is slightly smaller 
because of a thinner current sheet. This is still a first- 
order Fermi process, however. For comparison, we have 
also investigated the acceleration in a pure 3D turbu- 
lence, where the particles with gyroradii smaller than the 
injection scale accelerate through a second-order Fermi 
process. 

In summary, we have shown that the acceleration 
within reconnection sites, especially in the presence of 
turbulence, can be extremely efficient. This could be 
a powerful mechanism not only in the solar corona and 
wind or the Earth magnetotail, but also around black 
holes/accretion disks and jet launching regions of AGNs 
and GRBs. These results also call for further extensive 
work on the CR acceleration in magnetic reconnection 
sites with the inclusion of relevant loss mechanisms of 
CRs in order to assess the importance of this acceleration 
mechanism in comparison to other processes (e.g., diffu- 
sive shock acceleration) and to reproduce the observed 
light curves of the sources. 
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